A POLYNOMIAL GENERALIZATION OF THE EULER CHARACTERISTIC 

FOR ALGEBRAIC SETS. 

MIGUEL A. MARCO-BUZUNARIZ 

Abstract. We present a method to compute the Euler characteristic of an algebraic subset 
of C n . This method relies on clasical tools such as Grobner basis and primary decomposition. 
The existence of this method allows us to define a new invariant for such varieties. This 
invariant is related to the poblem of counting rational points over finite fields. 



1. Introduction 

One of the main invariants of a topological space is its Euler characteristic. It was initially 
defined for cell complexes, but several extensions have been defined to more general classes of 
spaces. In the setting of complex algebraic varieties, the natural extension is the Euler character- 
istic with compact support. In || Szafraniec gives a method to compute the Euler characteristic 
of a complex algebraic set by using methods from the real geometry. In this paper, we present 
another method, that only makes use of the basic properties of the Euler characteristic, and 
classical results on algebraic sets. This way of computing the Euler characteristic gives naturally 
a stronger invariant, which we define. 

The method works as follows. Consider V C C™ an irreducible algebraic set of dimension d 
and degree g. Take a generic linear projection it : C™ — > C d . If we consider tt restricted to V, it 
is a g : 1 branched cover. The branchig locus A and its preimage n \y X (A) can be computed. 
From the aditivity and the multiplicativity for covers of the Euler characeristic, we have the 
following formula: 

X (V) = g ■ x(C d ) - g ■ x(A) + x (tt I,; 1 (A)). 

So the computation of x(V) is reduced to the computation of the Euler characteristic of algebraic 
sets of lower dimension, allowing us to use a recursion process. 

In the previous method, we make use of the fact that x(C d ) = 1. If instead of making 
this substitution, we keep track of x(C d ) as a formal symbol, we obtain a stronger invariant 
F(V). This invariant is defined as a polynomial in Z[L], and has some interesting properties: 
the dimension, degree, and Euler characteristic of an algebraic set can be computed from this 
polynomial. It also gives information on the number of points on some varieties over finite fields. 
This relation with finite varieties could be used to compute this invariant by counting points. 
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In Sections |2| and |3| we show the preliminary results that prove the correctness of the method 
to compute the Euler characteristic, and describe the algorithm. Section |] is devoted to the gen- 
eralization of this method to the new invariant, which is also defined and some of its properties 
are shown. The extension of this invariant to projective varieties is discused in Section |^. In 
Sections^ and (7) we include an implementation of the two algorithms in Sage, together with some 
examples and timings. As an important example, we show that in the case of hyperplane ar- 
rangements this invariant coincides with the characteristic polynomial. Finally, the relationship 
of the invariant with the number of points over finite fields is shown in Section 

2. Theoretical justification 

Let V = V(I) C C™ be the algebraic set determined by a radical ideal /. Without loss of 
generality, we can assume that it is in general position (in the sense that we will precise later). 
By computing the associated primes of / we obtain the decomposition in irreducible components 
V = V\ U • • ■ U V c . The Euler characteristic x(^) can be expressed as x(Vi) — x({Vi) n (^2 U 
■ • • U V c )) + x{V 2 U ■ • • U V c ). The variety (Vi) H (V2 U • • • U V c ) is an algebraic set of lower 
dimension. So, by a double induction argument (over the dimension and over the number of 
irreducible components), we may reduce the problem of computing x(V) to the case where V is 
cither zero-dimensional or irreducible. 

If V is zero-dimensional, it consists on a number of isolated points, and its Euler characteristic 
equals the number of points. This number of points can be computed as the degree of the 
homogenization of the radical of / (which can be computed via the Hilbert polynomial, see |ffl, 
Chapter 5] for example). 

For the case of an irreducible variety V = V(I) C C™ being / < C[xi, . . . , x n ] a radical ideal 
of Krull dimension d 7 we will distinguish the homogeneous case from the non homogeneous. 

If I is a homogenous ideal, the variety V has a conic structure (it is formed by a union of lines 
that go through the origin). It means that V is contractible and hence its Euler characteristic 
is 1. 

For the non homogeneous case, consider the projection 

7T : C" ->■ C d 
(xi,...,x n ) n- (xi,...x d ) 

We may assume (aplying a generic linear change of coordinates if necessary) that the following 
condition is satisfied: 

Definition 1. Consider Ih < C[xo,xi, . . . ,x n ] the homogenization of /. We will say that / is 
in general position if y/I h + (x ,xi,. ..,x d ) 3 (£B , x\, . . . , x n ). 

Theorem 2. The previous condition is satisfied by any ideal I after a generic linear change of 
variables. Moreover, when this condition is satisfied the map it restricted to V is surjective and 
has no vertical asymptotes. 

Proof. If we consider the projectivization V C CP™, the projection n consists on taking as a 
center the n — d — 1 dimensional subspace S = {[xq : xi : ■ ■ ■ : x n ] | xq = X\ = • • • = Xd = 0}. 
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Since the dimension of V is d, the intersection V fl S is generically empty. We may hence assume 
that after a generic linear change of variables, S fl V = 0. This intersection is given precisely by 
the ideal y/Ih + (^0,^1, • • • ,x<i), which is homogeneous. The condition of V fl S being empty is 
equivalent to the ideal / being in general position. 

In this situation the preimage by tt |y of a point [1 : x\ : • • • : Xd] is given by the intersection 
of the subspace {[y : J/i : • • • : y n ] \ Vi = yoXi,---,Vd = yoXd} with V. By the genericity 
assumption, this intersection does not have points in the infinity. By dimension arguments, this 
intersection cannot be empty, and must be contained in the affine part of V. We have then 
proved that tt restricted to V is surjective. 

□ 



The intersection of a generic linear subspace of dimension n—d with V is a union of g distinct 
points, being g the degree of 1^. This degree can be computed through the Hilbert polynomial. 
Since / is in general position, all of the intersections of V with the fibres of tt will happen in the 
affine part. This means that tt restricted to V is a branched cover of degree g. We will see now 
that the branching locus of this cover is contained in a subvariety of C d that can be computed. 

Assume J = (/1, . . . , f s ) is in general position. Consider the matrix 



/ dfi 



M := 



dfi\ 



dfs 



df s 



and the ideal J generated by its (n — d) x (n — d) minors. 



Theorem 3. The branching locus of tv |y is contained in the elimination ideal (I 
C[ Xl ,...,x d }. 



J) n 



Proof. Consider a point p = (x\, . . . ,x n ) G V. If the linear space 7r _1 (7r(p)) intersects V at p 
transversally, then there is no ramification at p, since it means that at a neighbourhood of p the 
map 7T \ v is a difeomorphism. This condition of transversality can be expressed as follows: the 
normal space of V in p and the normal space of 7r _1 (7r(p)) generate the tangent space of C" in 
p. The normal space of V in p is generated by the rows of the matrix 



dfs 



(P)J 



The normal space of 7r _1 (7r(p)) is generated by the first d vectors of the canonical basis. A 
gaussian elimination argument tells us that these two spaces generate the whole space if and 
only if the matrix M has rank (n—d). So the set of points of V where tt |y ramifies is contained 
in the set S of zeros of / + J. 

The elimination ideal C[xi, . . . xa] (1 (I + J) is the Zariski closure of n(S). 

□ 
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Since V\tt 1 (tt(S)) is a cover over C d \ ir(S) of degree g, we have that 
X(V) = X(V \ TT-^niS))) + X (V n Tr-^TrfS))) = • X (C d \ tt(S)) + X (V D vr^MS))) = 

= 5 • (x(C d ) - x(tt(S))) + x(V n tt-H^SO)) = .9 - 5 • X (tt(S)) + X (V n tt-Htt^))). 

Both ir(S) and V PI tt^ 1 (tt(S)) are varieties of dimension smaller than V, so, by induction 
hypothesis, we can compute their Euler characteristic in the same form. 

3. Description of the algorithm 

Now we will describe, step by step, an algorithm to compute the Euler characteristic of the 
zero set of an ideal I = (fi, . . . , f s ). 

Algorithm 1. (Compute the Euler characteristic of the algebraic set defined by the ideal I): 

(1) Check if I is homogeneous. If it is, return 1. 

(2) Compute the associated primes (7i,...,J m ) of I. This can be achieved by primary 
decomposition (see Jl], Chapter 4]). 

(3) If there is more than one associated prime, we have that 

x{v(i)) = x(v(h)) + x(v{h n • • • n i m )) - x (v(h + {h n • • • n i m ))). 

by recursion, each summand can be computed with this algorithm. The following parts 
of this algorithm consider only the irreducible case, since we have already computed the 
associated primes, we will assume that I\ is prime. 

(4) Compute the dimension d and the degree g of V(I). If d is zero, return g. 

(5) Check that I is in general position. This can be done by computing a Grobner basis of 
y //, + (xo, xi, . . . , Xd) (where Ih is the homogenization of I) and using it to check that 
Xd+i, ■ ■ ■ , x n are in it. If it is not in general position, apply a generic linear change of 
variables and start again the algorithm. 

(6) Construct the ideal J generated by the (n — d) x (n — d) minors of the matrix 



/ dfi ... gfi \ 



M := 



i _dfs_ ... df 3 J 

\Xd+l Xn / 



(7) Compute the elimination ideal K = (I + J) (1 C[xi, • • • , xj\. 

(8) Compute by recursion x(V(K)) and x(V(I + K)). Return the number g — g ■ x(V(K)) + 
X(V{I + K)). 

4. A FINER INVARIANT 

The previous method essentially consists in decomposing our variety V in pieces, each of 
which is compared to C l through linear maps that are unbranched covers. At the end of the day, 
it gives us a linear combination (with integer coefficients), of the Euler characteristic of C\ 
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Now we will show that we can actually keep the information in this linear combination, defining 
a slightly different invariant. This information will be kept in a polynomial F W (V) £ Z[L], where 
L % will play the role of xi^ 1 ) 

We follow the same method as before but with two differences: 

• If the ideal I is homogeneous, we don't end returning a 1. Instead, we continue the 
algorithm, taking as Ih the ideal generated by I inside K[xq, . . . , x n ]. 

• In the final step, we return g ■ L d — g ■ F^{y{K)) + i^(y(7 + K))) instead of g — g ■ 
X(V(K))+ X (V(I + K)). 

So the algorithm results like this: 

Algorithm 2. (Compute the polynomial Fir(V) associated to an algebraic set V(I) in general 
position). 

(1) Compute the associated primes . . . , I m ) of /. 

(2) If there is more than one associated prime, we have that 

F n (V(l)) = F n {V{h)) + f w (v(i 2 n • • • n i m )) - F n [v{h + (j 2 n •• • n i m ))). 

by recursion, each summand can be computed with this algorithm. The following parts 
of this algorithm consider only the irreducible case, since we have already computed the 
associated primes, we will assume that I± is prime. 

(3) Compute the dimension d and the degree g of V(I). If d is zero, return g. 

(4) Construct the ideal J generated by the (n — d) X (n — d) minors of the matrix 

/ Ofi ... dfi \ 

Xd+l Xn 



M := 



V M*- ... dfjLl 



(5) Compute the elimination ideal K = (I + J) fl C[x l5 • • • , x^\. 

(6) Compute by recursion F n (V(K)) and F V (V(I + K)). Return the number gL d — g ■ 
F W (V(K)) + F„[V{I + K)). 

Note that both algorithms |l| and ^| can run differently if we apply a linear change of coordinates 
to / (which would change the projection ir). The topological properties of the Euler characteristic 
tells us that the final result of the algorithm [j] will coincide with the Euler characteristic regardless 
of this linear change of coordinates. But in the case of F^iV) we cannot ensure such a result. 
Nevertheless, for two sufficiently generic projections, algorithm ^ will follow the same exact steps, 
so we can define F(V) as the polynomial obtained by the algorithm]^ for generic projections. 

More preciselly, there must exist a Zariski open set T C GL(n, C) such that, the polynomial 
F 7r (a(V(I))) is the same for every linear change of coordinates a G T. 

Definition 4. Given an ideal I < C[x-y, . . . , x n ], we define the polynomial F(V(I)) as the 
polynomial F n (a(V(I)) for any a E T. 

We will say that I or V(I) are in generic position, or that we are in generic coordinates 
ifF w (V(I)) = F(V(I)). 



6 



MIGUEL A. MARCO-BUZUNARIZ 



Since so far we have no algorithmic criterion to determine if a projection is generic enough 
or not, the generic case can be computed by introducing the parameters of the projection, and 
computing the Grobner basis with those parameters. 

Anyways, experimental evidence suggests the following conjecture: 

Conjecture 1. If an ideal is in general position, it is also in generic position. 

Some partial results in this direction are easy to show. 

Lemma 5. If V in in general position, the leading term of F V (V) coincides with the leading 
term of F(V). 

Proof. It is immediate to check that the degree of ^(V) coincides with the dimension of V, and 
that the leading coefficient of F V (V) coincide with the degree of V, regardless of the projection 
used to compute it. □ 

Remark 6. The value of F n (V) at L = 1 equals x(V), independently of the choices of projections 
made for its computation, as long as we are in general position. 

These two results actually shows that F(V) is independent of the projection for the case of 
curves (since in this case it is a degree 1 polynomial whose leading term and value at 1 are fixed) . 
We will now show that the invariant F„ behaves well with respect to the product of varieties: 

Proposicion 7. Let I\ < C[xi, . . . , x n ] and I 2 < C[yi, . . . , y m ] be two ideals on polynomial rings 
with separated variables, and let V\ C C" and V 2 C C" 1 be their corresponding algebraic sets of 
dimensions d\ and d 2 respectively. Consider the ideal 

I ■= h + h < C[xi . . . ,x dl ,yi, ■ ■ ■ ,yd 2 ,x dl+ i, . . . ,x n ,y d2+ i, . . . , y m ]. 

Its corresponding algebraic set is V = Vi x V 2 C C" x C m = C n+m . Then F V (V) = F w (Vi) • 
F V (V 2 ). 

Proof. Without loss of generality, we can assume that we are in the irreducible case. We will 
work on induction over the dimension. If V\ or V 2 ar zero dimensional, the result is immediate. 

Consider gi,g 2 the degrees of V\ and V 2l and g the degree of V. It is easy to check that 
S = .9i ' .92- Is is also immediate to check that, if V 2 = C m , the statements holds (that is, 
F n (Vi x C m ) = F w (Vi) ■ L m ). Consider also Ai, A 2 and A analogously. 

Now we will show that A = (Ai x C d2 ) U (C dl x A 2 ). Let p = {x\, . . . ,x dl ,yi, ■ ■ -,y d2 ) & 
C dl x C d2 ). The set of points in V that project on p is the product of the set of points in V\ that 
project in {x\ . . . , x dl ) and the set of points in V 2 that project in (j/i, . . . , y d2 ). This set has less 
than gi ■ g 2 points if and only if (x\, . . . , x dl ) G Ai or (j/i, . . . , y dl ) e A 2 . It is immedaite also 
that (Ai x C d2 ) n (C dl x A 2 ) = Ai x A 2 . By induction hypothesis, we have that 

F W (A) = L d * ■ F w (Ai) + L dl • F^(A 2 ) - F n (A 1 ) ■ F V (A 2 ). 

Reasoning analogoulsy, we can conclude that 

^(tt-^A)) = F 7T (n- 1 (A 1 )) ■ F V (V 2 ) + F^n' 1 (A 2 )) ■ F^V,) - F^-^A^) ■ ^(^^A,)). 
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So, sumarizing, we have that 

FAV) = 3i 9 2(i dl+d2 -F 7r (A))+F 7r ( 7 r- 1 (A)) 
F V {V X ) = , 9l (L dl -F,(A 1 )) + F 7r (T r - 1 (A 1 )) 
F„(V 2 ) = g 2 (L d ^~FAA 2 ))+F n (TT~ 1 (A 2 )). 

Using all the previous formulas one can easily check that F^iy) = F n (Vi) ■ Fk(V 2 ). □ 

If conjecture [j] is true, the same result would be true for F(V). In fact, a weaker condition 
would be enough: if the product of two generc projections is generic, then the invariant F is 
multiplicative. This could be useful, for example, to give a criterion to check if an algebraic set 
can be the product of two nontrivial algebraic sets. If F(V) is irreducible in Z[L], then V couln't 
be a product. 

5. The projective case 

To compute the Euler characteristic of the projective variety V defined by a homogeneous ideal 
h <! C[#o, . . . , x n] we can also use algorithm [!} In order to do so, we will consider the hyperplane 
H "at infinity" given by the equation xo = 0. This allows us to decompose V as its affine part 
V := V \ H and its part at infinity V°° := V H H. It is clear that X (V) = x(V) + x(V°°)- 

The affine part V is an affine variety defined by the ideal obtained by substituting xo = 1 in 
the generators of Ih, whose Euler characteristic can be defined as seen before. 

The part at infinity V°° is a projective variety embedded in a projective space of less dimen- 
sion. The homogeneous ideal that defines is obtained by substituting X$ — in the generators 
of Ih- Its Euler characteristic can be computed by recursion. If we are in the case of CP 1 , V 
will consist on a finite number of points, which can be computed as the degree of \flh- 

Now we will show a different way to compute the Euler characteristic of a projective variety 
using the polynomial F(V). 

Theorem 8. Let I = C[xo, . . . ,x n ] ■ (/i,...,/ s ) be a homogeneous ideal in generic position. 
Assume that the generators /i, . . . , f s are homogeneous. Denote by Iq := (I + C[xq, . . . , x n ] ■ 
(xq)) fl C[xi, . . . , x n ], and l\ := (I + C[xq, • ■ • , x n ] ■ (xq — 1)) n C[xi, . . . , x n ]. That is, the ideals 
that represent the intersection ofV(I) with the hyperplanes {xq = 0} and {xq = 1} respectively, 
seeing the two hyperplanes as ambient spaces. Then the following formula holds: 

F(V(I)) = (L - 1) • F(V(h)) + F(V(I )) 

Proof. By induction on the dimension of V(I). If the dimension is zero, V(P) must consist only 
on the origin, since I is homogeneous. In this case, V(I\) is empty, and V(Iq) is also the origin. 
We have that 

1 = F(V(I)) = (L - 1) • + 1 = (L - 1) • F{V(h)) + F(V(I )). 

If the dimension d of V(I) is positive, consider the ideals J, K and H = I + K as before. 
Construct also Hq, H\,Kq and K\ in the same way as Iq and I\. Note that, since we are in generic 
position, the ideals H' and K' needed to compute F(V(Iq)) are precisely H and K (that is, 
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specializing xo — and then computing the minors of the matrix M and the elimination ideal 
is the same as computing the minors of the matrix and the elimination and then specializing). 
The same happens with J\ and K\. 
By induction hypothesis, we have that 

F(V(H)) = (L - 1) • F(V(H!)) + F(V(H )) 

and 

F(V(K)) = (L - 1) • F(V(Kx)) + F(V(K )). 

Now we have that 

F(V(I)) = g-L d -g-F(V{K))+F{V(H)) = 

= g ■ L d — g ■ ((L — 1) • F(V(Ki)) + F(V(K ))) + (L - 1) • F(V(Hi)) + F(V(H )) = 

= (L - 1) • (g ■ L^ 1 - g ■ F{V{K{)) + F(V(Hi)) + g ■ L^ 1 - g ■ F(K ) + F(H Q ) = 

= (L-l).F(V(h))+F(V(I )) 

and this proves the result. □ 

This theorem allows to relate the invariant F of the affine algebraic set defined by a homoge- 
nous ideal, and the invariant F of the projective variety defined by the same ideal as follows: 

Corolary 9. Let I be a homogeneous ideal in C[xi,...,x n ] in generic position. Let V be 
the algebraic set defined by L, and V the projective variety defined by the same ideal. Then 
F(V) = (L-1)-F(V') + 1. 

Proof. Aplying the previous result recursively, we have that F(V(I)) equals 

(L-iyF(V(h))+(L-l).F(V((Io)i))+- ■ .+(L-1).F(V((- ■ ■ (J„) • • • ) )i))+F(V((- ■ • (J ) • • • ) )). 

Note that if we compute piecewise -F(V') we obtain preciselly F(V(Ii)) + F(V((Io)i)) + ■ ■ ■ + 
F(V((- ■ ■ (Io) ■ ■ ■ )o)i))- Since V((- • • (L ) ■ ■ ■ )o)) consists only on the origin, we have the result. 

This last result can be interpreted as the fact that V is the complex cone over V. That is, 
V \ {0} is the product V x C* . 

□ 

6. Examples 

Both the polynomial F(V) and the Hilbert polynomial Pj have the same degree, and the 
leading term determined by the degree of V. This would point in the direction of considering 
that they contain the same information. The following example shows this is not the case: 

Example 10. Consider the conies Ci,C 2 G C 2 given by C\ := V(x 2 + y 2 - 1) and C 2 := 
V(x 2 + y 2 ). If we compute the Hilbert polynomial of the corresponding homogenous ideals in 
C[x, y, z] we get that P( x 2 +y 2 +z 2 ) = P^+yi) = 2 • t + 1. 

Let's now compute the polynomial V(C\) using the canonical projection tt : C 2 — > C to the 
first component. This projection is a 2 : 1 cover of C branched along the points ±1. Over each of 
these points, there is only one preimage. So finally we have that F{C\) = 2{L— 2) + 1 + 1 = 2L — 2. 
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On the other hand, the curve C2 also projects 2 : 1 over C, but now there is only one branching 
point (x = 0). So the result is F{C 2 ) = 2(L - 1) + 1 = 2L - 1. 

That is, this example shows that the polynomial F(V) contains information that is not con- 
tained in the Hilbert polynomial. 

An important example of algebraic sets is the case of hyperplane arrangements. We will now 
recall some related notions now (see [|[ Chapter II]). 

Definition 11. Let A be a hyperplane arrangement in C™. Its intersection lattice L(A) is the 
set of all its intersections ordered by reverse inclusion, with the convention that the intersection 
of the empty set is C" itself. 

The Mobius function is the only function // : L Z satisfying that: 

fi{C n ) = 1 

Ey<x»(Y) = via\{P}. 

Definition 12. The characteristic polynomial of A is defined as 

X (A,L):= ]T n{X)-L dim W. 

XeL(A) 

Theorem 13 (Deletion- Restriction). Let A be a hyperplane arrangement in <C n , and H a hy- 
perplane of A. Let A' be the arrangement resulting from eliminating H from A, and let A" 
be the hyperlane arrangement inside H induced by the intersection with A' . Then the following 
formula holds: 

X (A,L) = x(A',L)-x(A",L). 
Now we will see how this characteristic polynomial relates to the polynomial F (V). 
Theorem 14. Let A be a hyperplane arrangement in C™. The following holds: 

F(A) = L n - x (A,L). 

Proof. By induction on the number of hyperplanes. he case of only one hyperplane is immediate. 

If there are more than one hyperplane, take one hyperplane H of A. Since A — A' U H , we 
have, by aditivity, that 

F(A) = F(A') + F(H) - F(A' flfl) = F(A') + V 1 ' 1 - F(A"). 

Both A' and A" are hyperplane arrangements with less hyperplanes than A, so, by induction 
hypothesis the following formulas hold: 

F(A') = L»- X (A',L) 
F(A")=L n - 1 - X {A",L). 

Substituting these formulas in the previous one, and using the delition-restriction theorem we 
get the result. □ 

This esult tells us that the characteristic polynomial of A can be though of as the polynomial 
F of its complement. 



10 



MIGUEL A. MARCO-BUZUNARIZ 



7. Code and timings 
Here we show an implementation in Sage ([Q) of the two algorithms. 

7.1. Implementation of Algorithm [l]. 

def Euler_characteristic (I) : 
R=I.ring() 
if I . is_one () : 

return 
if R.ngens()==l : 

return sum( [j [0] . degree () for j in I . gen() . factor ()] ) 
J1=I. radical () 
if Jl . is_homogeneous () : 

return 1 
primdec=Jl . associated_primes () 
Jl=primdec [0] 
m=len(primdec) 
if m>l: 

J2=R. ideal (1) 

for j in [1 . .m-1] : 

J2=J2 . intersection(primdec [j] ) 

return Euler_charact eristic ( Jl) +Euler_char act eristic (32) -Euler_charact eristic ( J1+J2) 
P=J1 . homogenize () .hilbert_polynomial() 
if P . is_zero () : 

deg=0 
else : 

deg=P . leading_coef f icient *P . degree () . f actor ial () 
if deg==l : 

return 1 
dim=Jl . dimensionO 
n=R.ngens() 

varsl=R.gens() [0:n-dim] 
vars2=R.gens() [n-dim:n] 
varpiv=varsl [-1] 
IH=J1 . homogenize () 
S=IH.ring() 

JH=IH+S. ideal (S. gens () [n-dim:]) 
if JH. dimensionO >0: 
det=0 

while det==0: 

MH=random_matr ix (R . base_r ing ( ) , n) 
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det=MH . determinant ( ) 
L=list (MH*vector (list (R . gens ( ) ) ) ) 
return Euler_characteristic(R.hom(L) (Jl)) 
if dim==0: 

return deg 

M=matrix( [ [f . derivative (v) for v in varsl] for f in Jl.gensO]) 

J=R . ideal (M . minors (n-dim) ) 

K=(J+J1) . elimination_ideal(varsl) 

S=PolynomialRing(R. base_ring() , vars2) 

H=R.hom( [S(0) for j in varsl]+[S(j) for j in vars2] ) 

C=deg-deg*Euler_characteristic(H(K) )+Euler_charact eristic (K+Jl) 

return C 

This algorithm may be very slow (since it involves several Grobner basis computations), but 
in several interesting cases, it gives a useful answer in reasonable time. Let's show here some 
examples. 

In the case of curves and surfaces, the result is often reasonably fast; but it may vary a lot if 
a random change of variables has to be aplied. Here we show a few examples. These tests have 
been run on a Dual-Core AMD Opteron 8220. 

Three examples of plane curves: 

sage: R.<x,y>=QQ[] 

sage: time Euler_characteristic (R. ideal(x~5+D) 
5 

Time: CPU 0.16 s, Wall: 0.17 s 

sage: time Euler_characteristic (R. ideal(y~4+x~3-l) ) 
-5 

Time: CPU 0.19 s, Wall: 0.20 s 

time Euler_characteristic (R. ideal (x~2+y~2-5*x~2*y~4+x*y-l) ) 
-8 

Time: CPU 17.82 s, Wall: 17.82 s 

A curve and a surface in C 3 : 
S.<x,y,z>=QQ[] 

timeit ( ' Euler_characteristic (S . ideal (x~5+y~2+2*x*y+l , 3*x-5*y*x+y~2+l) ) ' ) 
10 

Time: CPU 0.17 s, Wall: 0.18 s 

timeit ( ' Euler_characteristic (S . ideal (x~5+y~2+2*x*y+l) ) ' ) 
-3 

Time: CPU 0.49 s, Wall: 0.49 s 

7.2. Implementation of Algorithm |^. 

@parallel(7) 
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@cached_f unction 
def FV(I,var='L') : 

FS=PolynomialRing (ZZ , var ) 

L=FS.gen() 

R=I.ring() 

if R.ngens()==0: 

return 
if I.is_zero() : 

return L~R.ngens() 
if I . is_one () : 

return 
if R.ngens()==l : 

return FS (sum( [j [0] . degree () for j in I .gen() . factor ()]) ) 
J1=I. radical () 
if Jl . is_homogeneous() : 

Sl=PolynomialRing(R.base_ring() ,R.gens() [0:-l] ) 

Hl=R.hom(list(Sl.gens())+[Sl(l)]) 

H2=R.hom(list(Sl.gens())+[Sl(0)]) 

resulp=FV([Hl(Jl) ,H2(J1)]) 

d=dict([[a[0] [0] [0] ,a[l]] for a in resulp] ) 

[il,i2]=[d[Hl(Jl)] ,d[H2(Jl)]] 

return (L-l)*(il)+i2 
primdec=Jl . associated_primes () 
Jl=primdec [0] 
m=len(primdec) 
if m>l: 

J2=R. ideal (1) 

for j in [1 . .m-1] : 

J2=J2 . intersection (primdec [j] ) 
resulp=FV( [Jl , J2 , J1+J2] ) 

d=dict([[a[0] [0] [0] ,a[l]] for a in resulp]) 
[il,i2,i3]=[d[Jl] ,d[J2] ,d[Jl+J2]] 
return il+i2-i3 
P=J1 .homogenize () .hilbert_polynomial() 
if P . is_zero () : 

deg=0 
else : 

deg=P . leading_coef f icient () *P . degree () . f actor ial () 
dim=Jl . dimensionO 
if deg==l: 
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return FS(L~dim) 
n=R.ngens() 

varsl=R. gens () [0:n-dim] 
var s2=R . gens ( ) [n-dim : n] 
varpiv=varsl [-1] 
IH=J1 . homogenize () 
S=IH.ring() 

JH=IH+S . ideal (S . gens ( ) [n-dim : ] ) 
if JH. dimension ()>0: 
det=0 

while det==0: 

MH=random_matr ix (R . base_r ing ( ) , n) 

det=MH. determinant () 
L=list (MH*vector (list (R . gens ( ) ) ) ) 
return FV(R.hom(L) ( Jl) ) 
if dim==0: 

return FS(deg) 

M=matrix( [ [f . derivative (v) for v in varsl] for f in Jl.gensO]) 

J=R . ideal (M . minors (n-dim) ) 

K=( J+Jl) . elimination_ideal (varsl) 

S=PolynomialRing (R . base_r ing ( ) , var s2) 

H=R.hom( [S(0) for j in varsl] + [S(j) for j in vars2] ) 

d=dict([[a[0] [0] [0] ,a[l]] for a in FV( [H(K) ,K+J1] )] ) 

[il,i2] = [d[H(K)] ,d[K+Jl]] 

C=deg*FS(L~dim-il)+i2 

return C 



This implementation makes use of the Sage framework for parallel computations, allowing to use 
several processor cores at the same time to compute the intermediate steps. It also caches the 
already computed results in case they would be needed later. 

Note that, if we assume Conjecture [l| to be true, this implementation works fine giving as entry 
the ideal whose algebraic set we want to compute (assuming it is in general position, which is 
something that can be easily checked). If we want to be safe from the posibility of the conjecture 
to be false, we have to introduce it over a ring that contains the parameters of the possible linear 
transformations. But, unluckily, the methods to compute the primary decomposition do not 
work over rings with parameters. One way to proceed is to compute the primary decomposition 
over the original ring without parameters. Then compute the discriminant with parameters, 
and then choose a value for the parameters in the open part of the Grobner cover (see g] for a 
definition and algorithm). 

Again, this method can be very slow, but in some cases it is fast enough to be useful. Here 
we present some of those examples: 
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The complex 2-sphere: 
sage: R. <x,y ,z>=QQ [] 

sage: time FV(R. ideal (x~2+y~2+z~2-l) ) 
2*L~2 - 2*L + 2 

Time: CPU 0.07 s, Wall: 0.32 s 

Another surface: 

sage: time FV(R. ideal (x~3+y~3+z~3-l) ) 

3*L~2 - 6*L + 12 

Time: CPU 0.07 s, Wall: 0.29 s 

The intersection of the two: 

sage: time FV(R. ideal (x~3+y~3+z~3-l ,x~2+y~2+z~2-l) ) 
6*L - 15 

Time: CPU 0.08 s, Wall: 43.01 s 

and their union (the timing is done after cleaning the cache of the function): 

sage: time FV(R. ideal ( (x~3+y~3+z~3-l) * (x~2+y~2+z~2-l) ) ) 

5*L~2 - 14*L + 29 

Time: CPU 0.08 s, Wall: 43.02 s 

notice the aditivity of the polynomial. 
The Whitney umbrella: 

sage: time FV(R. ideal (x*y~2-z~2) ) 
3*L~2 - 4*L + 2 

Time: CPU 0.13 s, Wall: 3.44 s 

The 3-sphere: 
sage: S . <x,y ,z , t>=QQ [] 

sage: time FV(S . ideal (x~2+y~2+z~2+t~2-l) ) 

2*L~3 - 2*L~2 + 2*L - 2 

Time: CPU 0.09 s, Wall: 0.43 s 

8. Counting points over finite fields 

In this section we will see how the polynomial P (V) can be related to the number of points 
of the variety considered over a finite field. Let's ilustrate this fact with an example. 

Example 15. Consider the conic given by the equation 

x^ ~\~ x^ 1 

in the affine plane over the field of 5 elements F5. 
The set of rational points is the following: 

(0,1), (0,4), (1,0), (4,0) 
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If we project it to the x\ axis, we see that over the point (0), we have two preimages, as 
expected by the degree. Over the points (1) and (4), we have just one point of the curve, since 
the cover ramifies there. But over the points (2) and (3) we have no points of the curve. The 
reason for this is that the equations 2 2 + x\ — 1 and 3 2 + x\ — 1 have no roots over F5. However, 
they do have all their solutions over a quadratic extension F25 of degree 2. In particular, if we 
look at the points of F5 x F25 that satisfy the equation, we obtain the set 

{(0, 1), (0, 4), (1,0), (2, a + 2), (2, 4a + 3), (3, a + 2), (3, 4a + 3), (4, 0)}, 

where a is an element of F25 that has minimal polynomial x 2 + 3 over F5. 

It might seem strange to consider the set of points in F5 x F25 that satisfy a given equation. 
But this set is in fact an algebraic set. Indeed, it can be expressed as the set of points of the 
affinc plane over F25 that satisfy the equations 

x^ ~\~ X2 1 
x\ — X\ 

Recall that the polynomial P(V) for such a conic was / = 2L — 2. In this case, we have 
obtained 8 points in total, which is preciselly the value of / for L = 5. A quick look at the points 
shows why this happens: there are two points over each value of F5, with the exception of the 
two branching points, where there is only one. 

Note that both algorithms |] and || can be run over finite fields in the same way as they 
run over the rationals. Whith a small exception, though: to ensure the existence of a change 
of coordinates that puts the ideal in general position, we might need to work in a finite field 
extension. Once done that, both algorithms would mimic the steps given by the algorithm run 
over Q, except if some leading coefficient becomes zero. But that would happen only for a finite 
number of prime numbers p. 

Notation 16. We can then define the polynomial F p (V) as the result of runing algorithm ^ in 
a field of characteristic p. 

As we have seen before, F p (V) — F(V) for almost every prime number p. 

Notation 17. Given a polynomial 

f = a + ai L + a 2 L 2 H h a n L n G Z[L] 

and a list of numbers (efi, . . . , d s ) with s > n, we will denote by f(d±, . . . , d s ) the number 

/(di, ...,d s )=a + aidi + a 2 did 2 H h a n did 2 ■ ■ ■ d n 

Another difference between the way the algorithms would run over Q and over finite fields 
lies in the primary decomposition. But then again, this will only happen for a finite number of 
primes. 
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Theorem 18. Given an ideal I < F p [xi, . . . , x n ], and the corresponding f := F P {V) of degree 
s, there exists a list of numbers (d\ < ■ ■ ■ < d n ) such that the number of points in 

(¥ pdl x ••• xF pJ „)ny(i) 

equals the number f(p dl , . . . ,p dn ). 

Moreover, for any such list (d\ < ■ ■ ■ < d n ) and any d\ multiple of di there exists another list 
(di, < ■ ■ ■ < di-i < d\ < d' i+1 < • • • < d' n ) satisfying the same property. 

Proof. By induction over the dimension. If / is zero dimensional, V(I) consists on F p (V) = 
deg(I) distinct points, whose coordinates lie in a suficiently big extension . Any sequence 
starting with d\ would satisfy the theorem. 

If I has dimension d > and degree g, we have that F p (V(I)) = g(L d — F p (V(K))) + 
F p (V(I + K)). By induction, we can assume that both K and I + K satisfy the result. Let 
(d\ < ■ ■ ■ < d\) and (d\ < ■ ■ ■ < d\) the corresponding sequences. Take d\ = gcd{d\,d\). There 
exist two sequences (d\ < d^' < • • • < d^) and (d\ < d%' < ■ ■ ■ < d^') that are valid for K and 
I + K respctivelly, and coinciding in the first term. Repeating this reasoning we can obtain two 
sequences (d\ < ■ ■ ■ < dd) and (d\ < ■ ■ ■ < da < d^+i <••• < d n ) that are valid for K and I + K 
respectivelly. 

Note that, since both F p (V{K)) and F p (V(I + K)) are of dimension at most d — 1, the terms 
dd, - ■ ■ ,d n can be changed arbitrarily and still the squence would be valid for K and I + K. 

Now take any point q := (qi,- ■ - ,qd) € (^ P d i x • • • x ¥ p d d ) \ V(K). Taking an apropiate field 
F p «g, we can ensure that there are exactly g points in {(xi, . . . ,x n ) £ (W p ™ q ) n \ x\ = qi, . . . , = 
Qd} H V. We can do the same for every point q, and take a common field extension F p « of all 
the different F p « g . This way, we have that there are exactly g points of V n (F^ x • • • x ¥ p d d x 
F p s x ■ ■ ■ x F p a) over each point of (F p<il x • • • x ¥ pdd ) \ V(K). 

By definition, we have that 

F p (V) = g(L d - F p (V(K))) + F p (V(I + K)). 

Now if we restrict ourselves to the points in ¥ p d 1 x • • • x ¥ p d d x ¥ pS x • • • x ¥ pS , we have that 

#V = g ■ (Z 1 / 2 ...p**- W(K)) + W(I + K). 

Making use of the induction hypothesis, and the fact that F p {V(K)) and F p (V(I + K)) are of 
dimension less than d, the result follows easily. 

□ 

Corolary 19. Let V be an algebraic set in C" defined by a an ideal I <J K[x\, . . . ,x n ], where 
K is an algebraic extension of Q. Then for almost every prime p there exists a list of positive 
integers (di iP , . . . , d ntP ) such that the number of points in 

(¥ pdl x • • • x F p d„) 

that satisfy the equations of I equals the number F(V(I))(p dl , . . . ,p d "). 
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This corollary could allow a different way to compute the polynomial F(V(I)) by counting 
points over finite fields. If we know the value of F(V(I))(S) for a sufficient number of such 
sequences S, recovering the coefficients of the polynomial is a simple linear algebra problem. 
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